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Photon orbital angular momentum and mass in a plasma 
vortex 
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Abstract. - We analyse the Anderson-Higgs mechanism of photon mass acquisition in a plasma 
and study the contribution to the mass from the orbital angular momentum acquired by a beam 
of photons when it crosses a spatially structured charge distribution. To this end we apply Proca- 
Maxwell equations in a static plasma with a particular spatial distribution of free charges, notably 
a plasma vortex, that is able to impose orbital angular momentum (OAM) onto light. In addition 
to the mass acquisition of the conventional Anderson-Higgs mechanism, we find that the photon 
acquires an additional mass from the OAM and that this mass reduces the Proca photon mass. 



Introduction. — Influenced by results derived in 1962 
by Schwinger [1], Anderson showed, in 1963, that a pho- 
ton propagating in a plasma acquires a mass, ^ 7 = huj p /c 2 , 
where u> p is the plasma frequency [2]. In this Anderson- 
Higgs process the photon acquires an effective mass be- 
cause of its interaction with plasmons, or better, a hidden 
gauge invariance in the plasma [3] . In order to study pho- 
tons that have acquired an effective mass, it is convenient 
to replace Maxwell's equations by Proca-Maxwell equa- 
tions [4,5]. In this Letter, we use this approach to analyse 
the contribution to the mass from the orbital angular mo- 
mentum acquired by a beam of photons as it traverses a 
spatially structured charge distribution. 

The orbital angular momentum (OAM) of light is 
a newly recognised observable of electromagnetic (EM) 
fields that is intimately related to optical vortices (OV's), 
phase defects embedded in certain particular light beams. 
OAM of light beams can be generated by the imprinting of 
vorticity onto the phase distribution of the original beam 
when it crosses inhomogeneous nonlinear optical systems 
[6] or particular spatial structures such as fork holograms 
or spiral phase plates. Such beams can be mathematically 
described by a superposition of Laguerre-Gaussian (L-G) 
modes characterized by the two integer-valued indices I 
and p, or by Kummer beams [7]. The azimuthal index I 
describes the number of twists of the helical wavefront and 



the radial index p gives the number of radial nodes of the 
mode. 

The EM field amplitude of a generic L-G mode, in a 
plane perpendicular to the direction of propagation, is 



F P i(r,(p) = 



Awpl 



obeying the orthogonality condition 
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where w is the beam waist and L™ is the associated La- 
guerre polynomial. The phase factor exp(— Hip) is associ- 
ated with an OAM of IH per photon, and a phase singular- 
ity is embedded in the wavefront, along the propagation 
axis, with a topological charge I [8,9]. The intensity dis- 
tribution of an L-G mode with p = has an annular shape 
with a central dark hole where the intensity vanishes be- 
cause of total destructive interference. 

As is well known, not only the linear momentum of light 
but also its angular momentum can propagate to infinity 
[10-12]. The OAM property of the field remains stable 
during the propagation in free space and has been exper- 
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imentally verified down to the single photon limit [13]. It 
has also been studied theoretically [14]. 

A photon in vacuum has its intrinsic angular momen- 
tum, the spin (SAM), that can assume only two values, 
u = ±1 corresponding to transverse left-hand and right- 
hand circular polarisation, respectively. If the photon were 
massive, then it would carry a third independent longitu- 
dinal polarisation component along the direction of propa- 
gation [15]. The exchange of angular momentum between 
a photon beam and a plasma vortex and the possible ex- 
citation of photon angular momentum states in a plasma 
was analysed in Ref. [16]. There can also be an exchange 
of angular momentum between electromagnetic and elec- 
trostatic waves in a plasma due to stimulated Raman and 
Brillouin scattering processes [17]. The properties of plas- 
mons carrying OAM were analysed in Ref. [18]. 

In this Letter we show that the OAM acquired by a pho- 
ton in a spatially structured plasma can be interpreted 
as an additional mass-like term that appears in Proca 
equations. More specifically, we study the propagation 
of a photon with wavelength A in a static helicoidally dis- 
tributed plasma with step qo = X/p where p is an integer 
[17]. This apparent mass term shows that structured spa- 
tial and temporal inhomogeneities of matter distribution 
can impose properties onto single quanta. This can be 
interpreted as a manifestation of a Mach principle in clas- 
sical and quantum electrodynamics. Henceforth in this 
Letter, we use natural units e = c = h= G = 1. 

Proca equations and OAM in a plasma vortex. 

The Maxwell-Proca Lagrangian density A describes the 
kinetic potential of a massive EM field in which there ap- 
pears a mass term p 2 A fl A fl /2 because of the chosen gauge 
invariance: 



currents p and j, are 
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(3) 



where /j" 1 is the reduced Compton wavelength associ- 
ated with the photon rest mass and j M = (p, — j) is the 
4-current. F^ v is the electromagnetic tensor and A^ the 
4-vector potential. From the Lagrangian density, one ob- 
tains the covariant form of the Proca-Maxwell equations 



7^ + ^ = ^ 

that leads to the Proca wave equation for A^ 
(□ - ntyAp = -47rj„ 



(4) 



(5) 



When one considers a photon propagating in a plasma, 
the usual Maxwell equations can be replaced by the set of 
Proca-Maxwell equations in which there appears a mass- 
like term for the photon due to light-matter interaction 
[4] . The usual formulation of Proca-Maxwell equations is 
obtained by expanding Eq. (4) in terms of the electric E 
and magnetic B fields, that, in the presence of charges and 



V • E = Airp — p 2 ^ 
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V-B = 

V x B 



(6) 
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and for p 1 — > 0, they smoothly reduce to Maxwell's equa- 
tions. The Poynting vector for massive photons depends 
directly on both the scalar and vector potentials 



-(ExB + u^A) 



(7) 



and also the energy density has an explicit dependency on 
the potentials 



u=±- (E 2 
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While the Lorentz invariance remains valid, the gauge 
invariance is lost because the potentials become observ- 
able through the acquired energy densities /x 2 <^ 2 /8/T and 
p 2 A 2 /87r, respectively, due to the interaction of the pho- 
ton with the plasma or its confinement within the observ- 
able Universe. 

Let us consider a transverse EM wave propagating 
through an isotropic plasma, cast to form a hclicoidal 
static plasma vortex. The heavy ions constitute a neu- 
tralising background, and their motion can, in the first 
approximation, be neglected. The motion of free electrons 
forms a three-dimensional current j = — nv, where n is 
the electron number density, and v the velocity of the 
electrons in the medium derived from the electron fluid 
equations 



dn 

~dt 
dv 

~dt 



V • nv = 



+ v 



Vv = (E + v X B) 

m 



(9) 
(10) 



Thermal and relativistic mass effects are ignored. In the 
first approximation, the mean electron velocity in the 
static plasma vortex becomes v = v (r, t) +6v where v is 
the background velocity and 5v is the perturbation associ- 
ated with the propagating EM wave. The electron number 
density is 



n = n a + n(r, z) cos(l ip + q z) 



(11) 



where uq is the background plasma density and the he- 
licoidal density perturbation in the plasma is given by 
the latter term, expressed in cylindrical coordinates, r = 
(r, <p, z). The coordinate z is the axis of symmetry around 
which the electron spiral is winding and along which the 
EM wave is propagating. 

Similarly as in a spiral phase plate, the number of elec- 
trons affecting the EM wave depends on z and can vary 
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slowly on a scale much larger than the spatial period 
zq = 2n/qo, where qo is the helix step. For a typical dou- 
ble vortex one obtains Iq = 1. Neglecting for a moment 
the rotation of plasma and considering the case of a static 
helical perturbation, the current density perturbation in 
the plasma is j = — n(r)5v, and the propagation equation 
of the electric held becomes 



dt 2 



w- IE 







(12) 



where uj 2 — w 2 [l + e(r, <p, z)) and cj p0 is the unperturbed 
plasma frequency, given by ui 2 — Aim^/m. The quantity 
e represents the perturbation. We consider solutions of 
the form 



E(r, t) = a(r) cxp ( - iuit + i J k(z') dz'^j 



(13) 



where uj is the EM wave frequency and a(r) is the ampli- 
tude, varying slowly along z such that [17] 



d 2 a 



dz 2 
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(14) 



We can then write the wave equation in a perturbed parax- 
ial equation form 

(Vi + 2ife^-w=)a = (15) 

and the dispersion relation, connecting k and to, is 

k 2 = io 2 -u 2 p0 [l + e(r,<p,z)} . (16) 

In our case, a general solution to the wave equation in 
the paraxial approximation can be represented in a basis 
of orthogonal L-G modes, having the amplitude 

a(r, <p, z) = b P i (r, z)e d * exp ( - ^ ) e pl (17) 
P i 

where w = w(z) is the beam waist, e p i are unit polarisa- 
tion vectors, A p i are the amplitudes of each mode and 



where the function l) p represents the associated Laguerre 
polynomial. As usual, the integers p and I are the radial 
and the azimuthal quantum numbers, respectively. The 
electric field then becomes 



E 

with 



(r,t) =^E pi (r)cxp(-kjt + i f fc(z')d/) (19) 
P i J 



E p( (r) = c p i(z)F p i(r,(p)e p 
where F p i(r, tp) is given by formula (1). 



(20) 



When a vortex perturbation e(r, ip, z) is present, the 
modes will couple according to 

§- z c A z ) = ^J2 K ^p' 1 '^' ( 21 ) 



where K(pl,p'l') are the coupling coefficients. If the EM 
wave does not carry OAM and the mode coupling is suffi- 
ciently weak that the zero OAM mode dominates over the 
entire interaction region, one obtains a coupling 

K(pl,p'l') = ™ p 2 — V[^,- io e i9oZ +Si>,i e- i «>*] (22) 

If we assume the same polarisation state for all the inter- 
acting modes, the field mode amplitudes describe the rate 
of transfer of OAM from the static plasma vortex to the 
EM field 



c P,±kXz) 



7tc(0) r u 2 (z') n(z') 
2c 2 J k(z') n 



e^ iqoZ ' dz' . (23) 



A more general solution, where the amplitude of the ini- 
tially excited mode is allowed to change, is discussed in 
Rcf. [16]. The initial OAM state of the electromagnetic 
beam decays into other states (/, + ulo) on a length scale 
approximately determined by the inverse of the coupling 
constant, showing an effective exchange of OAM states 
between the photons and the plasma. 

The photon mass in Proca equations is defined as to 7 = 
fiy and the effective photon mass is also related to the 
plasma frequency m e s = oj p [2]. By comparing these two 
definitions, one obtains the equivalence between the in- 
verse of the characteristic length in a plasma, /i 7 and 
the plasma frequency, = uj p . When considering the 
EM wave equation of the Proca field one obtains a Klein- 
Gordon equation for the 4-vector potential 



(□ - ul)A^ = -Airj^ 



(24) 



with the constraint derived from the massive photon in a 
medium (d^A^ = 0). By differentiating this expression 
with respect to time, and considering, in addition to the 
previous calculations, the simplest case where /i 7 is a con- 
stant in time, one obtains 



(□ 



d . 



A 9 ■ 



(25) 



We now consider the spatial components applied to the 
case of photons moving in a plasma. In particular, we 
consider a plasma with a well-defined structure of a static 
plasma vortex, in which the density n = n(r) is not a 
function of time. Being the set of equations independent 
on the temporal part, we pass from the four-dimensional 
notation to a three-dimensional vectorial notation without 
losing in generality. The current density has the form 
j = — n(r)v, and its derivative becomes 



d_ . 

dt ]l 



^v^n(r) + n(r)^v) = - n (r)^v (26) 
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and the wave equation of the Proca EM field becomes 
d 



(□-/4)^A = -47m(r)^v. 



d_ 

di 



(27) 



Using E = - 
the EM field 



— Mr A one obtains the wave equation for 



(□-/£)(e + v« 



AnnM—v. 

K ' dt 



(28) 



Assuming that v is parallel to E, and v = v + <5v, where 
5v is the perturbation of electrons velocity associated with 
the propagating EM wave, we obtain for E ^ the wave 
equation for the EM field in the plasma, 



v • V<5 



-47T 



n(r)Sv — v • D"V(f) 

m 



E = (29) 



where v = v/\v\ is the direction vector of the velocity field, 
v the unit vector of velocity and Sv = v-dtv = \dtv\. When 
comparing this equation with Eq . (12), derived from the 
electric field propagation equation in the case of a static 
plasma vortex, one obtains 



M7 



1 + 



E 



+ 47T 



n(r)Sv 



□ Vc 



E 



(30) 



that implies a direct relationship between the effective 
mass that a photon acquires in a plasma, the plasma fre- 
quency and the orbital angular momentum because of the 
peculiar spatial distribution: 



a4 = 



E 



E + v-Vcf,^ 

- -= r = , ( 4n5v \n + n cos(l ip + q z 

E + v • V0 V 

- 4ttv • nVcj^j 



)] (31) 



The new mass component is a fictitious term, generated 
by the interaction of photons with the plasma and cannot 
be ascribed to an intrinsic property of the photon. 

When the electron number density exhibits certain spa- 
tial properties, such as vortices, any photon has an associ- 
ated virtual mass term that is smaller than that expected 
from Proca equations in a homogeneous plasma, because 
of a negative term that corresponds to a precise orbital 
angular momentum component. 

Conclusions. — We have investigated the problem of 
photon mass in a plasma and show that part of the ac- 
quired mass term is related to the orbital angular momen- 
tum of light imposed by certain spatial distributions of the 
plasma. We focused our attention on the simplest case of 
spatial distributions described by a static plasma vortex. 
This approach shows that also the spatial distribution of 
charges can impose OAM and an additional mass term 
that reduces the effective mass of the photon inside a non- 
structured plasma, degrading the variation in OAM. The 



spiral-like plasma structure induces OAM states in most of 
the scattered photons and the stochastic interaction value 
responsible for the Proca mass term is instead transformed 
into an organised state of light, with the result of reducing 
the averaged mass term [19]. 
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